Theory of spin noise in nanowires 
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We develop a theory of spin noise in semiconductor nanowires considered as prospective elements 
for spintronics. In these structures spin-orbit coupling can be realized as a random function of 
coordinate correlated on the spatial scale of the order of 10 nm. By analyzing different regimes of 
electron transport and spin dynamics, we demonstrate that the spin relaxation can be very slow 
and the resulting noise power spectrum increases algebraically as frequency goes to zero. This 
effect makes spin effects in nanowires best suitable for studies by rapidly developing spin-noise 
spectroscopy. 
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Nanostructures are the promising hardware elements for 
■spintronics [l[ - a rapidly developing branch of physics 
and technology aiming at studies and application of spin- 
dependent phenomena in the charge transport and infor- 
mation processing. The quest for the systems with ul- 
■tralong spin relaxation times d is one of the main chal- 
lenges in this field. Since the dynamical spin fluctuations 
Q characterized by correlations on the spin relaxation 
timescale, are seen as a spin noise in the frequency do- 
main, this search can be done with recently developed 
highly accurate low-frequency spin noise spectroscopy ji| 
aimed at the measurement of intrinsic equilibrium spin 
dynamics. The spin noise spectroscopy allows to study 
•the slow spin dynamics in (llO)-grown quantum wells [5[ 
and in quantum dots |(|. Theoretical background of this 
method is given, e.g., in Refs. 

An interesting class of semiconductor nanostructures 
demonstrating peculiar and slow spin dynamics are the 
quantum wires [Ml, where e.g. InAs, InSb as well as 
GaAs/AlGaAs systems are the prospective realizations. 




The effects of spin-orbit (SO) coupling on the transport 



pup 

were clearly demonstrated there [ 1 3 Jl4( and the nanowire 



.based qubits were introduced jl5[|lo|. A SO coupling in 
duced effective magnetic field acting on electron spins in 
nanowires is directed parallel or antiparallel to a certain 
'axis [l7 - 21 1 resulting in a giant spin relaxation anisotropy 
similar to that expected in some two-dimensional sys- 
tems [22| . Since the SO coupling is a structure- and 
•material-dependent property, all sorts of disorder (ran- 
dom doping interface fluctuations [27[, random 
variations in the shape, etc.) which cause electron scat- 
tering and nonzero resistivity, can cause local variations 
in the coupling. As a result, in addition to the regular SO 
coupling, caused by the lack of bulk (Dresselhaus term) 
or structure (Rashba term) inversion symmetry, all low- 
dimensional structures inevitably have the random con- 
tribution in it . The spatial scale of the fluctuations is of 
the order of 10 nm as determined by the characteristic 
distances in nanostructure were shown to give rise to a 
number of fascinating phenomena 2^-33- However, their 
role in quantum wires was not studied so far. 

Here we address theoretically the electron spin dynam- 
ics in ballistic and diffusive semiconductor nanowires aim- 



Figure 1: Schematic plot of the experimental configuration: a 
quantum wire (dark stripe) is illuminated by a linearly polar- 
ized beam and Kerr rotation angle of its polarization plane 8k 
is measured. Polarizations of the beams are marked by double- 
headed arrows. Dashed arrow corresponds to the polarization 
of the reflected beam in the absence of the Kerr effect. 



ing at the study of the spin noise spectrum. Different 
regimes of electron spin relaxation are determined and 
the crossovers between them are analyzed in detail. In 
particular, we demonstrate that when the electron mo- 
tion is diffusive and the dominant contribution to the SO 
interaction is random, the spin relaxation becomes alge- 
braic rather than exponential and the spin noise power 
spectrum diverges at low frequencies to as 1/w 1 / 2 , show- 
ing colored noise [31 - 33| well suited for the studies by the 
spin noise spectroscopy. A very slow spin dynamics result- 
ing in the low-frequency noise divergence makes nanowires 
an exception among semiconductor systems. 

The spin noise spectroscopy, reviewed in Ref. is 
based on the optical monitoring of the spin fluctuations 
[33 | in Faraday, Kerr or ellipticity signals measured with 
a weak linearly polarized probe beam incident on a single 
wire or a wire array sample, see Fig. Q] It can be shown 
similarly to Refs. [J, 0, [35| that for the probe tuned to 
the fundamental absorption edge, the Kerr rotation angle 
Ok oc s z , hence its autocorrelation function is directly 
related to the spin noise: (9K(t)0K{t')) oc (s z (t)s z (t')) , 
where s z (t) is the density of the z— component of the total 
electron spin. As a result, this optical technique measures 
long-range correlations of equilibrium spin fluctuations oc- 
curring in the illuminated spot. 



We consider a single channel quantum wire extended 



2 



along the a;— axis and represent the SO Hamiltonian as: 



T-Lso = -[a(x)k x + k x a{x)]a x . 



(1) 



Here k x = —id/dx is the electron wave vector component 
along the wire axis, a(x) is the coordinate-dependent SO 
coupling strength. In Eq. ([T]) we assumed that the spin 
quantization axis, A, is fixed, and a\ is the component of 
spin operator along this axis. The specific form of the SO 
Hamiltonian Eq. (JlJ implies that the effective field acting 
on electron spin points either parallel or antiparallel to 
the axis A. This is obvious for a constant a{x) (l7-21|. 
and holds true provided that the microscopic symmetry 
of the fluctuations forming the SO coupling randomness 
is the same as overall symmetry of the system. 

The SO coupling is assumed to be the sum of the 
coordinate- independent contribution, ao, and the Gaus- 
sian random function with zero average, a r (x) such as 
a(x) = ao + a T (x) with the correlation function [29| : 



(a r (x)a r (x')) = (a^)F co „(x - x'), 



(2) 



where (pq) is the mean square of SO coupling fluctuations 
and the range function F corr (x — x'). We introduce also 
the typical correlation length of the SO coupling 



In 



F corr (x)dx, 



(3) 



characterizing the size of the correlated domain of the 
random SO coupling. Details of the models of random 
SO coupling can be found in Ref . (29l| . 

We begin with the semiclassical regime, where SO cou- 
pling disorder is smooth on the scale of electron wave- 
length, Id ^> Ap, where the wavelength of the Fermi level 
electrons Af = 27r/fcp, with kp being the Fermi wave vec- 
tor for the degenerate electron gas. The Hamiltonian ([TJ 
implies that the spin rotation angle around the A-axis dur- 
ing the motion from the point xq to x\ is 



2m f Xl 
9(xi,x ) = — a(x')dx', 

h Jxn 



(4) 



where m is the electron effective mass. Eq. (j4|) shows that 
the angle is solely determined by electron initial and final 
positions and does not depend on the history of the motion 
between these points. This result, being well established 
for the systems with regular SO coupling [2ol [36l - l39l | holds 
also for the nanowires with the SO coupling disorder. As 
it follows from Eq. ([I} the spin precession rate is pro- 
portional to the electron velocity and given coordinate- 
dependent function. Hence, it does not matter whether 
the electron starting from the point xq reached the point 
X\ ballistically or diffusively: all contributions to spin pre- 
cession of the closed paths, where electron passes the same 
configuration of a(x) in the opposite directions, cancel 
each other. 

The temporal evolution of electron spin is directly re- 
lated the electron motion along the wire. We consider 
here spin projections at given z axis, perpendicular to the 
spin quantization axis A. Time dependence of electron 



spin z component averaged over its random spatial mo- 
tion and over the random precession caused by the field 
a(x) can be most conveniently characterized by the cor- 
relator (s z (t)s z (0)) = (sl(0))C ss (t) with the normalized 
correlation function: 



C ss (t)= dxp(x,£)(cos[6»(x,0)]), 



(5) 



where p(x, t) is the probability that electron travels dis- 
tance x during the time t. Note, that C ss (t) can be un- 
derstood as disorder-averaged electron spin z component 
found with the initial condition s z (0) = 1. It results from 
the linearity of the spin dynamics equations: the correla- 
tors (sj(i)sj(O)) satisfy exactly the same equations as av- 
erage values (si(t)} = x,y,z). In derivation of Eq. 
we assumed also that the scattering of electrons, which 
determines p(x, t) is not correlated with the random SO 
field oii(x), hence, the averaging over the realizations of 
a r (x) denoted by the angular brackets and over the tra- 
jectories can be considered independently. This can occur 
in nanowires where random Rashba fields are induced by 
doping while the momentum scattering is due to the wire 
width fluctuations. If the same local disorder determines 
the electron scattering and random SO fields, in relatively 
clean systems the electron mean free path I exceeds by far 
the disorder correlation length Id in Eq.Q. Hence, spatial 
scales of two random processes: I for the electron back- 
ward scattering in the random potential and Id for the 
spin precession are strongly different. As a result, on the 
Z-scale, the memory of the short-range correlations is lost, 
and Eq. ([5]) holds. Although Eq. ([5]) is presented for the 
smooth SO coupling disorder, where the electron motion 
is semiclassical, Z^/Af 3> 1, a general Green's function 
approach confirms it for arbitrary Id/ Ap values. 

Our next step is to perform averaging of cos [6{x, 0)] in 
Eq. ([5]) over the random realizations of the a(x)-field. For 
this purpose we recast 



cos [0(x,0)} 



( * — fj?^~ x 



exp [h9 r (x) 



where 



tf r (x) = 2m/H 2 / a r (x')dx' 



(6) 



(7) 



is the contribution of the random SO coupling into the 
spin rotation angle. We expand last exponent in series 
in $ r assuming the Gaussian SO coupling disorder. In 
the averaging, odd powers of spin rotation angle vanish, 
(0 2n+1 (x)) = 0, for integer n and even powers can be 
expressed solely with (6^(x)) as 



<<r (*)> 



2m 



a r (x')dx' J = (2n-I)!!(6» r 2 (x))". 

(8) 

Direct calculation shows that the mean square (6f(x)) 
caused by the random SO interaction is given by 



9?(x))=2^) 2 (a r 2 )£dx'£ dyF a 



r (y). (9) 



3 



Finally, Eq. ([5]) reduces to 



C„(t) 



A I 4-\ ( 2ma ° 

ax p(x, t) cos — —5 — x 



cx P [-(^(.t))/2]. 

(10) 

When (9 2 (x)) becomes considerably larger than one, spins 
are completely dephased. Equation (fTU|) is our central re- 
sult: it relates temporal average spin dynamics with elec- 
tron motion along the wire. Distribution function of elec- 
tron displacements, p(x, t), presented for different regimes 
of electron motion below, enables us to calculate spin evo- 
lution by Eq. ([T0| . The spin noise power spectrum is given 
by the transform of C ss (t) Q : 



(s 2 z ) u = 2 C 3S (t) cos (uit)dt. 



(11) 



To get a better insight into the problem, we begin with 
the key limits [4l[. First, for the ballistic electron dynam- 
ics p(x,t) = S(x — VFt), where vf = HIcf/th is the Fermi 
velocity. The ballistic motion is realized on the temporal 
scale t <C t = l/vp with r being the momentum relax- 
ation time. We are interested in the spin dynamics on 
the time scale t ^> Td = Id/vp, where T<j is the time dur- 
ing which electron passes the correlated interval of the SO 
coupling fluctuations. Using Eq. ([TO]) we obtain damped 
oscillations of the spin z-component: 



C ss (t) sa cos (Q t) exp (-t/r s>r ), 



(12) 



with the frequency Qq = 2mctoVF/h 2 determined by the 
averaged SO coupling and the decay time caused by the 
SO coupling fluctuations 



1 



2mvF 



K 2 )r d . 



(13) 



Equation (|13p for the spin relaxation time t st is a result of 
random spin precession (29[. Spin noise power spectrum 
calculated using Eqs. (fTT|) and (fT12j) reads: 



(4\ 



2t s r Re 



1 



fi 2 T 2 
12 T s, 



r + (1 - iujT S ^y 



(14) 



with the result presented in Fig. [2j 

This ballistic regime of spin dynamics, however, can be 
realized only in very clean systems, where f2oT 3> 1. Oth- 
erwise, electron spin evolution occurs at the time scale, 
where electron moves diffusively (Fig(31 upper panel) , i.e. 



p{x,t) = 



1 



-x 2 /4Dt 



2^fnDt 



(15) 



where D = v f t is the diffusion coefficient. In the 
absence of the SO coupling fluctuations and provided 
that Qqt *C 1 exponential spin relaxation is due to the 
Dyakonov-Perel' mechanism 17], |2l| with the relaxation 
time t S: dp = l/(^o r )- The spin noise spectrum has a 
Lorentzian form (s 2 .)^ = 2t Si dp/(1 + w 2 t s 2 dp ) with the 
width determined by the relaxation time. 

New physical features arise when the SO coupling fluc- 
tuations dominate over the regular contribution. From 




Figure 2: Spin noise power spectrum, (s z )u, for ballistic prop- 
agation, fioT"s,r = 2. Due the exponential decay in Eq. (|12|) it 
is finite at tJ = with the width determined by the spin re- 
laxation time T Sir . The spectrum peaks at the frequency fio 
since average electron spin rotates in the SO field at the rate 
fio and asymptotically decays as lj~ 2 in accordance with the 
fluctuation-dissipation theorem. 



now on we put ao = and consider the system where SO 
coupling is purely random and concentrate on the long- 
time (t ^> Td 7 T 7 T S:I ) dynamics. At these times, the sys- 
tem in Eq. (|10[) is characterized by two length parameters. 
One parameter is the diffusion length y/Dt in Eq. (fl~5j) , the 
other one 



dxcxp [-(e 2 {x))/2\, 



(16) 



characterizes spin randomization. At sufficiently long 
times, when \[~Dt ^> L s , one can take p(0,t) instead of 
p(x, t) and immediately obtain from Eq. (|10|) that the re- 
laxation is algebraic rather than exponential: 



C„(()«p(0,i) / da; exp [-{9 2 {x))/2\ = 



(17) 



Equation (|17[) predicts extremely long spin dccohcrence 
described by the inverse square root law: (s z (t)) cx \j\ft. 
This surprising result has a transparent physical inter- 
pretation (see Fig. [3]) : Indeed, if an electron is displaced 
from its initial position by a sufficiently large distance, 
x > L s , its spin rotation angle becomes so large, that 
it does not contribute to the total spin polarization ow- 
ing to exp [-(9 2 (x)) /2] in Eq. (fTU|) . As a result, the spin 
polarization is supported by the electrons located in the 
vicinity of their initial positions, mainly due to the return 
after multiple scatterings by the random potential. The 
fraction of such electrons, in agreement with the diffu- 
sion distribution, decays as p(0, t) oc resulting in the 
same behavior in the spin polarization. It is interesting 
to mention that this qualitative argument does not work 
for the constant SO coupling despite spin of electron is 
restored upon the return to the origin also here. The rea- 
son is that due the oscillations of the spin on the spatial 
scale of the order of ti 2 /mcto (see Fig. G2 lower panel) in 
Eq. fTO)) , the diffusive return of electrons to the origin is 
insufficient for formation of the algebraic relaxation tail. 

Another realization of the \j\fi spin decay can be 
achieved for the very strong random SO couping where 
the spin relaxation occurs within one nanosize domain of 
the SO coupling, that is at the electron displacement much 
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Figure 3: Upper panel: Schematic illustration of the dis- 
placements distribution p(x, t) for two different time moments: 
to < t\. Lower panels: The quantum wire and spins of diffus- 
ing electron for the random and regular SO couplings, respec- 
tively. For the random SO coupling, if the electron is within 
the L 3 distance from its initial point [see Eq. I|17|l] its spin is 
preserved, when it leaves this interval, the spin dephases. 



less than Id- In this case, spin relaxation rate is due to 
the Dyakonov-Perel' mechanism and is determined by the 
local value of a(x) inside the domain. Spins of electrons lo- 
cated in the intervals with large a(x) will relax fast, while 
spins of those experiencing weak a(x) will relax slow. 

Slow non-exponential spin relaxation, described by 
Eq. (|17p manifests itself in the low frequency spin noise 
spectrum. From Eq. (fTTj) it follows: (s^) oc i.e 
the spin noise diverges at oj — > 0. Such a non-trivial be- 
havior is inherent to the quantum wires with random SO 
coupling, where spin restores upon return to the origin: in 
multichannel wires for sufficiently fast intcrchannel scat- 
tering (42l | and in two-dimensional systems spin relaxation 
is exponential [29[ and ( s i) u - is finite. 

To conclude, we studied theoretically spin noise in semi- 
conductor nanowirc for different regimes of the electron 
propagation. We demonstrated that if the spin relaxation 
is determined by the randomness in the SO coupling, spin 
relaxation becomes algebraic being closely related to the 
high probability for electron to stay close to its initial po- 
sition as a result of a multiple scatterings in the random 
potential. This behavior can appear in at least two possi- 
ble regimes: (i) when the electron motion is diffusive and 
(ii) when the spin relaxation occurs on a small spatial 
scale of the order of 10 nm. In any of these cases, the spin 
noise power spectrum shows colored l/y/uj noise. In ad- 
dition, this observation shows that low-frequency optical 
spin noise spectroscopy is an excellent tool for studying 
spin phenomena in semiconductor nanowircs and charac- 
terization of random potential and SO coupling there. 
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Supplementary Material for "Th eory of Spin Noise in Nanowires" 



SI. SPIN NOISE AT ARBITRARY FREQUENCIES 

Here we determine the spin noise spectrum for arbitrary 
frequencies w. We employ the kinetic equation for electron 
distribution function f(x,v x ,t) dependent on the position, 
velocity v Xl and time. The equation has the form: 



dt 



df | /-/ 
dx T 



0, 



(SI) 



where f(x, v x , t) satisfies the initial condition f{x, v x , 0) = 
d(x)[S v ^ VF + S Vt - VF ]/2 meaning that at t = it is built at 
x = with the equal fractions of electrons with veloc- 
ities v x = ±vp. As a result, the carriers can be sepa- 
rated into the right movers, v x = vp, and left movers, 
v x = — v-p with function / = [f(x,vp,t) — f{x, — vp, t)]/2 
being the anisotropic part of the distribution. The dis- 
tribution of electron displacements is given by p{x, t) = 
f(x, vp,t) + f(x, —vp, t). It can be shown that the spatial 
Fourier transform and cositot) transform of this distribu- 
tion, p(k,uj) has the form 



p(k,u) = 2 Re 



t(1 — iujt) 



(kl) 2 — \lot{1 — \wt) 



(S2) 



In accordance with Eq. (11) in the main text the spin 
noise spectrum can be presented as 



(* 2 X = 



dk 
2tt 



p(k,cj)T(k), 



where 



/oo 
exp [ikx- (6 2 {x))/2]dx. 
-OO 



(S3) 



(S4) 



Analytical result can be obtained in the regime where 
spin rotation angles within each correlated domain of the 
SO coupling are small, that is VL?Td = 2m^/ (a 2 )ld/h <C 1 
with f2 r = 2yJ (a 2 )k x /h. Here the spin dynamics occurs 
on the spatial scale i» Id, mean squares of spin rotation 
angles are proportional to electron displacement (0 2 (x)) ps 
2(£l r Td) 2 \x\/ld being valid for x 3> Id or at t 3> Td, and 
function T(k) takes the form: 



T(k) 



2l d (n r T d ) 2 
(n r r d ) 4 + {ki d y 



(S5) 



After lengthy transformations we obtain 

/ 2\ f.-^ T d 



2Rc 



(n r Td) 2 y/iu)T/(iuJT - 1) - lUJTd 



(S6) 



It can be seen from Eq. (S6) that at low frequencies, u <C 



St" 1 



spin noise spectrum has the form: 



(4\ 



V2T.T 



(S7) 



in agreement with the analysis above. For high frequencies 
lot ^> 1, {s 2 )^ is given by 2/{uj 2 t s ^) since at r d <C t <C r 
the electron motion is ballistic, and spin dephasing is 
caused by the random fluctuations of the spin-orbit cou- 
pling, cf. Eq. (12) of the main text. The entire frequency 
dependence of (s^) is plotted in Fig. |4] 



P 

i 

3 

A 4 

Nn 10 



CO 

V 




10" 4 0.001 0.01 0.1 1 
Frequency (dt) 



Figure 4: Spin noise power spectrum for diffusive electron 
propagation, Qo = 0, Q Y Td = 0.01, rj/r = 0.1. Solid line 
shows exact result, calculated according to Eq. (S6). Dotted 
(with the slope -1/2) and dashed (with the slope -2) lines show 
the low-frequency and high-frequency asymptotic, respectively. 



SII. SPIN DYNAMICS AND NOISE IN 
MULTICHANNEL WIRES WITH RANDOM 
SPIN-ORBIT COUPLING 

The spin evolution in multichannel structures depends 
on the additional set of parameters, {Tjj} being the scat- 
tering times between the channels i and j, as well as on the 



G 



details of spin dynamics in every channel. For qualitative 
analysis (a general case requires a separate treatment) we 
consider a structure with two conducting channels, where 
(i) spin-orbit coupling disorder in different channels is not 
correlated and (ii) in each channel Q r Td <C 1 (super- 
scripts denote channels) , i.e. spin rotation angles in corre- 
lated domains of the spin-orbit coupling are always small. 
Here we can characterize the interchannel scattering by a 
single time r c and focus on the most interesting case with 
no regular contribution to the spin-orbit held: ao = 0. 

In the limit of very rare interchannel scattering events 
(the condition is given below) the channels are indepen- 
dent. Hence, the general results expressed by Eqs. (10), 
(11) of the main text and by Eq. (S6) as well as asymp- 
totic Eqs. (S7) and (17) of the main text hold. Although in 
these equations one has to average over the realizations of 
a[ c '(x) in different channels, the low-frequency spin noise 
power spectrum remains \/y/uj, the same as in a single 
channel wire. 

Now we turn to the efficient interchannel scattering with 
short r c . If t c <C Td electron quits given channel faster 
than it quits the correlated domain. Spin rotations be- 
tween interchannel scattering events are uncorrelated and, 
due to this randomness, spin dynamics is exponential: 



(s z (t)s z (0)) ocexp(-r c t), 



(S8) 



where the 



relaxation 

2 



rate T r is of the order of 



max 



.. Similar exponential decay of the 

spin correlator remains for < t c < t. Here, the mean 
square of the spin rotation angle between interchannel 
scatterings can be estimated as ((<S<I>) 2 ) = (9^(v^t c )) cx 
t c , with v>- c ' being the characteristic velocity in the chan- 
nel. Spin relaxation is governed by the Dyakonov-Perel'- 
like mechanism, with the rate 



roc mf)/T c 



(S9) 



which is r c -indepcndcnt for exactly the same reason as the 



spin relaxation rate due to the random spin-orbit [Eq.(13) 
in the main text] does not explicitly depend on the elec- 
tron free path. 

Most interesting physics appears for a very weak in- 
terchannel scattering, r c ^> r, Td- In this case, electron 
moves diffusively in a given channel before the interchan- 
nel scattering occurs. As we have shown above, the 1/y/t 
tail in the spin polarization (and corresponding 1/y/tH spin 
noise) results from the carriers dwelling around the initial 
point of their trajectories. Since the tail is formed at long 
times \f~Di ^> L s , sec Eq. (16) of the main text, it is sup- 
ported by electrons which moved many times back and 
forth in the random potential. If the interchannel scat- 
tering is probable, electron may return to the initial point 
via other channels, where its spin rotation is not correlated 
with that in the initial one. Therefore, in general 1/Vt 
tail is destroyed and the usual exponential spin relaxation 
takes place. However, if r c is long enough to assure that 
the typical electron displacement during the diffusion be- 
tween interchannel scatterings L c = \fDr c 3> L s , there 
is a time interval L 2 S / D <C t <C r c and the corresponding 
frequency range, where the spin dynamics and the noise 
are algebraic: 



(*,(t)«,(0)} cx -L ( s 2) w cx _L (S10) 



In the regimes of a highly efficient interchannel scattering, 
with the spin relaxation described by Eqs.(S8) or (S9), the 
probability of spin components restoration upon return to 
the initial position is strongly suppressed, and, as a result, 
spin noise power spectrum at low frequency decreases and 
becomes finite. 

For completeness, we mention that if the random spin- 
orbit coupling docs not depend on the channel, that is 
Qi (x) = a[ 2 ' (x) , spin precession angle between any points 
Xi and X2 is insensitive to the interchannel scattering, and 
our analysis in the main text holds exactly the same. 



